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Primordial magnetic fields from second-order cosmological perturbations: 

Tight coupling approximation 
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We explore the possibility of generating large-scale magnetic fields from second-order cosmological 
perturbations during the pre-recombination era. The key process for this is Thomson scattering 
between the photons and the charged particles within the cosmic plasma. To tame the multi- 
component interacting fluid system, we employ the tight coupling approximation. It is shown 
that the source term for the magnetic field is given by the vorticity, which signals the intrinsically 
second-order quantities, and the product of the first order perturbations. The vorticity itself is 
sourced by the product of the first-order quantities in the vorticity evolution equation. The magnetic 
fields generated by this process are estimated to be ~ 10"'^^ Gauss on the horizon scale at the 
recombination epoch. Although our rough estimate suggests that the current generation mechanism 
can work even on smaller scales, more careful investigation is needed to make clear whether it indeed 
works in a wide range of spatial scales. 
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INTRODUCTION 



Magnetic fields are known to be present on various 
scales in the universe [1]. For example, magnetic fields 
are observed in galaxies and clusters, with intensity 
~ 1 /iGauss. Only an upper limit has been given for 
magnetic fields on cosmological scales, < 10~^ Gauss. 
Primordial large-scale magnetic fields may be present 
and serve as seeds for the magnetic fields in galaxies and 
clusters, which are amplified through the dynamo mech- 
anism after galaxy formation. We need numerical cal- 
culations in order to evaluate precisely a necessary seed 
magnetic field, however we can estimate order of magni- 
tude of it and it is ~ lO^^^ - lO^^^ Gauss at lOOMpc 
[2]. If sufficiently large, the primordial magnetic field 
could leave observable imprints on the cosmic microwave 
background [3]. 

A number of models have been proposed for generat- 
ing large-scale magnetic fields in the early universe [4- 
6]. However, they rely more or less on some unknown 
physics. In the present paper, we discuss magnetogene- 
sis in the pre-recombination era using only the conven- 
tional physics that has been established. The genera- 
tion of magnetic fields in this era has been studied in 
Refs. [7-15]. Now it is widely accepted that large-scale 
cosmological perturbations, generated from inflation in 
the early universe, evolve into a variety of structures such 
as the cosmic microwave background anisotropies and 
galaxies. However, inflation produces only density fluc- 
tuations (scalar perturbations) and gravitational waves 
(tensor perturbations); vector perturbations, and hence 



large-scale magnetic fields, are not generated in the con- 
text of usual infiationary scenarios, unless, for instance, 
modifying Maxwell theory of electromagnetism. Even 
if they were generated due to some mechanism, they 
only decays without any sources. This no-go argument 
is based on linear perturbation theory, and so we will be 
studying second-order perturbations to overcome this dif- 
ficulty within infiation and standard Maxwell theory. We 
consider a multi-fiuid system composed of photons, elec- 
trons, and protons, which are tightly coupled via Thom- 
son and Coulomb scattering but slightly deviate from 
each other [16]. Thomson scattering is important for 
the generation of large-scale magnetic fields in the pre- 
recombination era, because a rotational current will be 
produced by the momentum transfer due to the Thom- 
son interaction. We shall see how this process occurs by 
doing the tight coupling expansion. 

The time scales of Thomson scattering. Coulomb scat- 
tering, and the plasma oscillation at z ~ 10^ are tt = 



lO'^ sec. 



Tc = 10~^ sec, and u!~^ = W^" sec, respec- 
tively [14]. The last two time scales are much smaller 
than the Thomson time scale, implying that the effects of 
Coulomb scattering and plasma oscillations are smoothed 
out within the Thomson time scale. For this reason, we 
will not take seriously Coulomb scattering and plasma 
oscillations [11, 14]. And the effect of the synthesis of 
hydrogen, i.e., the ionization rate, on magnetogenesis 
may be taken into account carefully [15]. It has been 
reported that the change of the ionization rate does not 
affect the generation of magnetic fields on cosmological 
scales, though the effect can no longer be ignored on 
scales smaller than AU. Thus, for our current purpose, it 
is safe to ignore the effect of the synthesis of hydrogen. 

This paper is most closely related to the work of 
Ref. [11], which studied how magnetic fields are gener- 
ated from nonlinear cosmic inhomogeneities using the 
tight coupling approximation. While the covariant for- 
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malism was used in Rcf. [11], wc use in this paper cos- 
mological perturbation theory and so the two works are 
complementary to each other. Recently, Takahashi et al. 
took a simplified Newtonian approach to see the detailed 
physical processes that result in magnetogenesis in the 
pre-recombination era [14]. The present study may be 
viewed as the general relativistic extension of Ref. [14]. 
Second-order vector perturbations have also been stud- 
ied previously in different contexts [17 20]. Although 
our rough estimate suggests that the current generation 
mechanism can work, more careful investigation is needed 
to make clear whether it indeed works in a wide range of 
spatial scales. 

The organization of this paper is as follows. In the next 

section we define our key quantities and derive the equa- 
tions of motion based on second-order cosmological per- 
turbation theory. In Sec. 3, wc invoke the tight coupling 
approximation to solve the relevant equations, giving a 
compact formula describing the generation of magnetic 
fields. In Sec. 4 we give a summary of the present paper. 



II. BASIC EQUATIONS IN SECOND-ORDER 
COSMOLOGICAL PERTURBATION THEORY 

A. The energy-momentum equations up to 
second-order in perturbation theory 

We start with defining the metric and the energy- 
momentum tensor we use. Throughout the paper we 
rely on the energy-momentum (non)conservation equa- 
tions and so we explicitly write them down in this sec- 
tion. 

The background spacetime is given by the spatially 
flat Friedmann-Lemaitre-Robertson- Walker metric. We 
write the perturbed metric in the Poisson gauge as 

= a'^(ri) [- (1 + 2(6) drj"^ + 2xidr)dx' 

+ {1- 2n) Sijdx'dx^] , (1) 

where a is the scale factor and r] is the conformal time. 
We study the cosmological perturbations up to second 
order, and hence we write 

<^ = </,w + </.(2), xi = xf\ n = n('^+n^'\ (2) 

where we have dropped the first order vector perturba- 
tions since they are not generated from inflation in 
the standard scenarios. We also neglect the tensor per- 
turbations (gravitational waves) for simplicity. 

We consider a multi-fluid system composed of pho- 
tons (7), electrons (e), and protons (p). We assume that 
the energy-momentum tensor for each fluid component is 
given by that of a perfect fluid (i.e., we neglect anisotropic 
stresses): 



where pi = wipi (wp = = 0,Wy = 1/3) and is 
the 4-velocity of the fluid satisfying g^''u(^i)^U(^i)^ = —1: 
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'+ o'^^+ o^(/b-^(7)'^(/) ) ' 



2'^ -2^WjV)' 



V{i)i + Xi- 27^^;(/)i) . 
The energy-momentum tensor is explicitly written as 



rp 

J (7)0 



-pi - {pi +Pl)V(i)jV\ 



(4) 

(5) 

(6) 
(7) 



= {PI+PI){V^^I)^+X^) 

-{pi+Pi){cl) + 2n)v^i)i, (8) 



and 



'"{1)1^ = {Pi +Pl)V(I)^v'^j-^ +PiS/ 



(9) 



Note here that the energy density and pressure of each 
fluid, pj and p/, already include inhomogeneous parts: 
Pi = Pi{t- x) and pj = pi{t. x). We will do the perturba- 
tive expansion of these quantities later. 

The energy-momentum tensor of the electromagnetic 
field is written in terms of the field strength F^i, as 



-'EM 



Xcr 



(10) 



As usual, the field strength is decomposed into the elec- 
tric and magnetic fields as 



(7)1/) 



(11) 
g. Here we 



where e^"^ = 7]^""^" U(^^)„ with rp^"^^ 
tacitly adopt the photon frame. 

The vorticity, lo\ = —if 
portant quantity in our discussion. The vorticity is given 
explicitly by [9] 



ie"'''^W(/)AV,,?i(/)p, is the im- 



■■>i 

,(2)» 



1 \^(I) 



(12) 



(2)\ 

Xk') 



+<-^^K^+27e«)], (13) 

where Cijk is the 3-dimensional flat alternating tensor. 
Since we are assuming that the first order vector pertur- 
bations vanish, we have Wj^'* = 0. 

During the pre-recombination era photons and charged 
particles are tightly coupled via Thomson scattering, and 
charged particles are tightly coupled to each other via 
Coulomb scattering. The momentum transfer between 
fluids I and J are specified by sf ^. The momentum 
transfer due to Thomson scattering is given by 



hi)^ = iPi 



+ Pi)u(i)^,u1i)+pi5^ (J = 7,p,e), (3) 



-^arrieP^ (m(7)j - M(e)i) , 

ml At \ 



(14) 
(15) 



3 



where (Tt, P7, ^^e, np and TOe, TOp are the Thomson cross 
section, the energy density of photons, the number den- 
sity and mass of electrons and protons [11]. Each number 
density is defined in the rest frame. The difference be- 
tween frames in the momentum transfer terms appears 
at third order (in the form of the Lorentz factor), and 
hence can be neglected. The momentum transfer rate fol- 
lows from the collision term in the Boltzmann equation. 
A careful and detailed derivation of Eqs. (14) and (15) 
(valid at second order) can be found in Refs. [14, 21, 22]. 



The momentum transfer due to Coulomb scattering is 

written as 

= -nf = -e^npHeVc {u(p)i - U(e)i) , (16) 

where rjc is a electric resistivity which comes from 
Coulomb scattering [14]. To derive the energy conser- 
vation law and equations of motion, we compute the di- 
vergence of the energy momentum tensor: 



v.r(/)o" = - [p'l + 3Hipi + pi)] + 3{pi + pi)n' - [pi + pi)d,vlj^, (17) 

V.T(j)i" = [(P/ +P/)(W(/)^ + X^)]' + m{pi+Pl){V(i)i + Xi) + {pi+Pl)di{4) - </>2) + d^pi 

- (<^ + 2n) { [{pi +pi) V(i)i]' + m {pi+pi) «(/)/} - Hpi + Pi)'R'v(i), + ip,: + Pi)d, (^'(lypi^ (18) 

I 



where the prime denotes the derivative with respect to 
r] and H = a' /a is the conformal Hubble rate. Note 
that the energy transfer by Thomson scattering may be 
ignored [23] [8, 11]. Thus, the equations of motion gov- 
erning the present 3-fluid system are given by 



= k7 



eUpEi + Hi- 



pe 



-erieEi + K 



ep 



(19) 
(20) 
(21) 



where we neglected the Lorentz force from the magnetic 
field because it will give rise to higher order contributions. 
To proceed, it is convenient to define the center of mass 



velocity and the relative velocity as 



V{b)i ■ = 



dv, 



(22) 
(23) 



where (3 := rrie/mp. Hereafter we will assume charge 
neutrality: rip = ne(=: n) in the rest frame of charged 
particles. This will be justified because the velocity dif- 
ference between protons and electrons, (5w(pe)i, is ignored 
compared to Sv(^yb)i in the situation which the Coulomb 
interaction is very strong. 

Now the equation of motion for the photon fluid reads 



1^ P'-yi'^{i)i + Xi) + Pi{v{i)t + Xi)' + 4Hp7(v(7)i + Xi) + \dtPi -{4> + 27?.) Piv[^)i + P'iV{^)i + W.p^vi^^)i 



\aaTnp~f{l — 2TZ) 



The equations of motion for the charged particles (20) and (21) are combined to give 

+ P) { n'{i\h)i + Xi) + n'{v^b)i + Xi)' + 4:T-Cn{v(^b)i + Xi) - {(p + 27^) nv[^^^ + n'v(jj)i + AUnv^i^) 



(1 + /32) 5v(^b), + l^5v^p,)l2A) 



— bnR-'v, 



+ndi {(j) - (jp) + (w(p)iW(p) + /3u(e)it^(e)) > = jaarnp^il - 211) 



(1 + 0^) Sv^^b)i + T;^^^^(pe)i 



(6)i 



• (25) 



The equations of motion for the charged particles can be combined in a different way. The manipulation niene x (20) 
-nipUpX (21) leads to another independent equation: 

mempnj {nSv^pe)i)' + ^Un5v(^pe)i -{(t> + 2TZ) {nSv^pe)i)' + 4:nn5v(pe)i - 5n'n6v^pe)i + ndj (t^(p)jt^(p) - '?^(e)i^^(e)) | 
= en'^mp (1 + f3)Ei- ^aarn'^mpp^il - 27^) [(l - (3'^) 5v(^^t,)^ + '^^(pe)*] " e^mpU^l + (3)r]c{l - 27^)5^;(pe)i(26) 
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Since the Coulomb interaction is very strong in the 
regime we consider, we may take the hmit ^ 
Sv(pe)i — > 0. Then Eq. (26) simply reduces to 



1 - Aut 
1 + /3 3e 



ap-y(l — 2TZ)Sv, 



(27) 



This may be regarded as "Ohm's law" in some sense. In 
the standard Ohm's law, the electric field is proportional 
to the electric current density ~ eJiij^pgjj. It is reminded 
that the contribution from the electric current gives us 
the diffusion term in the evolution equation for the mag- 
netic field, and then the source for the magnetic field 
cannot be induced. However, the electric field is propor- 
tional to Sv(^jb)i in the above formula. The current in our 
"Ohm's law" is originated from the velocity difference 



between protons and electrons through the interaction 
with photons. Indeed, if one takes the same mass limit 
oi P ^ 1 (though it is not realized in the nature), the 
electric field cannot be generated. 



B. Magnetic fields 

Let us turn to the evolution of the electromagnetic 
fields. The Bianchi identities, V[A-F)ti/] = 0, yield the 
induction equation 

{a^B')' = -e'i'^dj [a (1 + <!>) Ek] - e'^^ {avjEk)' . (28) 
Using "Ohm's law" (27), this can be written as 



1-/33 Aar 
+ /? 3e 



I 



(29) 



It is now clear from Eq. (29) that we need the velocity locity difference can be computed as follows. The manip- 
difference between the photons and the charged particles, ulation (4p^/3)~^ x (24) — [nip (1 + /?) n]~^ x (25) leads 
to evaluate the generated magnetic field. The ve- to 



— +Xi) - — {'Vib)i + X^) + {Sv{yb}i)' + AU5v(^i,)i -{(t> + 2TZ) —V(^)i - —V(^b)i + {Sv{jb)iy + AH5v(^b) 



Pi 



n 



(76) 



+ 



1 diPy 
4 p^ 



Pi 



+ dj (^^(7)i^^(7)) - Y^^J (^(P)»^(p) + /5^(e)i^(e)) = ""(^ ' ^T^)SV(^ib)i, (30) 

r 



where a is defined as 



III. TIGHT COUPLING APPROXIMATION 



1 + 



with 



and 



3mr, 



3(1 + /3) 



R:-- 



3mp (1 + /3) n 

4/07 



TT 



aaTPj 



TT , 

(31) 



(32) 



(33) 



In deriving the above equation we again used ^ 
dv(jh)i- Thus, our remaining task is to solve Eq. (30) for 



dv 



"(7&) 
(7b)i- 



We are to solve Eq. (30) using the tight coupling ap- 
proximation (TCA) [16]. In this approximation the time 
scale of Thomson scattering {tt = rrie/ auTP-y) is as- 
sumed to be much smaller than the wavelengths of the 
perturbations (fc~^). Thus, the small expansion param- 
eter of the TCA is fcrr, which is dimensionless. During 
the pre-recombination era, photons, protons, and elec- 
trons are strongly coupled via Thomson scattering, and 
hence the TCA will be a good approximation. 

At zcroth order in the TCA, all fluid components have 
the same velocity Vi and the density fluctuations are adi- 
abatic. Following Ref. [14], we define the deviation from 
the adiabatic distribution for baryons by 



nb =nb{l + Afe) . 



Then, we expand various quantities such as A;, and 



(34) 
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in terms of the tight couphng parameter kTr- 



*6 



6v 



(76)1 



6v^ 



(b)i 

(I) 
(76) » 



{II) 



(11) 



+ • 



(35) 
(36) 

(37) 
(38) 



where Vi is the common velocity of photons and baryons 
in the tight couphng hmit. Our notation is that Roman 
and Arabic numerals stand for the order of TCA and that 
in cosmological perturbation theory, respectively. Here, 
we adopt the photon frame, so that A^^^ = A^^^' = 0. In 
other words, the quantities associated with photons give 
the "backgroimd" in the TCA. Note that we consider 
cosmological perturbation theory and the TCA simul- 
taneously. The following analysis includes cosmological 
perturbations up to second-order and the tight coupling 
expansion up to TCA(II), where TCA(n) denotes the 
tight coupling approximation at n-th order. The nota- 
tion X^"^") indicates that X is a m-th order quantity in 
cosmological perturbation theory at TCA(n). 

At each order in the TCA (n = /, //), Eq. (29) reduces 

to 



(») 



1 + 13 3e 



(n) 



(39) 



where 



(n) 



ijk 



(0) \°Pj °^{lb)k ) 

1 

+a,(<^(i)-27e(^))5.;-;) 



"3'^^(^b)k + _(0)' 

Pi 



+- 



-(0) V Pi ^3 '^^ilb)kj 



(40) 



We shall now evaluate the source term 

A. TCA(O) 

1. Energy conservation 
At (0,0) order, Eq.(17) reads 
(pf)'+4Hpf =0, (n(°))'+3W)=0. (41) 

At next order in cosmological perturbation theory, i.e., 
at (0, 1) order, we similarly have 

-(5p«)'-4H5p« + 4pW(7^(i)y 

-^pWay(i)=o, (42) 

-(^nW)'-3H<5nW + Sn^ (7^«)' 



These equations can be put into 



n(0) 
implying that 



^ = 3 (n^^^y - dy^^\ (44) 



3 V 



(45) 



2. Vorticity 

We show that the vorticity vanishes at TCA(O). Using 
Eq. (41), the momentum conservations for photons and 
baryons at (0,1) order are 



4 ) 



(46) 



We eliminate {v^^)' from the above equations and obtain 
the relation 



(1) _ ldi5^ 



(47) 



At (0,2) order, we eliminate similarly {vl +Xi )'• The 
momentum conservation gives 



( (2) , (2) 



l)+27^«')^;(^) 



1 



Taking the curl of the above equation and using Eq. (47), 
we can show that w^^^' = 0. It indicates that the vorticity 
is always zero at TCA(O). 



B. TCA(I) 

1. Energy conservation 

To compute Sv'^!yl^^\ we evaluate A^^'^^ using the en- 
ergy conservation. At (/, 1) order, Eq. (17) reduces to 

- (n^A^^'i))' - 3WnWA(^'^) - n^°^d,v[lf' = 0, (49) 
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and hence 



or, equivalently, 



which is integrated to give 



(50) 



(51) 



vP = - 



did pip 



+ ^(0)"' 4(1 + ^(0)) ■ 

(56) 

Taking the curl of Eq. (53), we arrive at the vorticity 
conservation equation (for the photon fluid): 



2. Vorticity evolution equation 

To derive the vorticity evolution equation, we take the 
curl of the total momentum conservation. The total mo- 
mentum conservation is 



(2a2pfa;'(2)y + 8a2Wpfa;'(2)=o 



(57) 



+PTdi{<j) -</>') + ^ P/(l + wi)dj (v^i)ivlj^^ = 0, (52) 



where we used Eq. (56). It is important to note here that 
there is no source for the vorticity at TCA(I). Therefore, 

(2) 

we may assume that ujy at TCA(I). 

3. Derivation of Sv'^^^^'^l and <5v|:^J,^'. 

At (/, 1) order, Eq. (30) reads 



^(0) n(0) 



(1) , 1 diSpj _(Q) / I 

» 4 -(0) (76)1 ^ 

P7 



where pb = mp{l + j3)n. Note that the scattering terms 
are cancelled out. At TCA(I), this equation is simplified 
to 

\pT{Vi+Xi)^ + ^'HpT{Vi+Xi)+PTdi{<t>-(j)'^) 

+ \diP^ -{<i> + 27e) [{pTVi)' + AHpTVi] 



and so we have 



{■yb)i 



1 



- ^prll'vi + prdj (viV^) = 0, 
where pt is defined by 

4 4 
Pt ■■= ^p-y + Pb= 3^7(1 + -K)- 

At (0, 1) order, Eq. (53) reads 



(53) 



(54) 



(1) 



1 d.Spi^^ 



(59) 



where 



a 



(0) _ (l+/3^)(l + i?("))4aaTpf 



1 + 



3m^ 



(60) 



^ Pt 



3 



= 0, (55) 



At next order in cosmological perturbation theory, i.e., 
at (/, 2) order, we have 



7 



l(0) 



n(0) n(0) 



p-f 



n' 



(0) 



„(i) 



+ 



where 



(P"f)' (n(0))' 



Pt^ 



7,(0) 



pt^ 



(61) 



v(l) 



1 



(l + ^(0)) 

(I 2) 

This can be solved for (V^J,)^ to give 

,(^.2) _ 1 



Pt^ 



5v 



v(0) 



v(i) 



me)Sn 



(1) 



4 + 3^W _(o)<5py^ 
4(1 + ^(0))" ' 



I P7 



-(0) 

Pt 



v(0) 



-(0) 

pV 



2 Pt^ 



(62) 



(63) 



C. TCA(II) 



yielding 



1. Derivation of Sv'^^P^'.^^^ and 5v\' ' '^"^ 



Eq. (30) at order (/, 1) and (/, 2) allows us to obtain 
the expression for (5w|;^^jj'' and respectively. At 

first order in cosmological perturbation theory, we have 



(n(o))' 



n' 



(0) 



(76)1 

(64) 



Sv 



(7b)i 



1 



(65) 



At second order in cosmological perturbation theory, we 
get 



J 



From this we obtain 

6)i - (7«')i^ ^ "''^(76)iJ + (a(0))2 [ (,''^(7fe)i j + ""^C 



"^(76)* 



"(76)i 



a(0) ^^(76)i 



+ - 



1 



(67) 



where we used Eq. (65) and note that 



2. Vorticity evolution equation 



a 



a 



(0) 



1 4- P(0) '~^b ~ 1 _i_ p(0) / '^''1 ^i^'^ljb) •(^^) In order to derive the vorticity evolution equation at 

i + H + n J TCA(II), we need to consider Eq. (52) at TCA(I). At 



8 



(/, 1) order, Eq. (52) reads 
, (1),, HR^°^ (1) 



1 



-a 



1 + ^(0)"' 4(l + i?(0)) p(0) 



1 + ^(0) 



(76)) 



where wc used Eqs. (41) and (54). We take the curl of 
Eq. (52) at (/, 2) order, so that we obtain 



(a2a;(2)')' + 



0) 



1 



2(l + ii(0)) 
where we used Eq. (69) and 



e^^'=a,a(^'^)<5^;|;,^)i (70) 



rv(0) 



(71) 



Since the right hand side in Eq. (70) contain the source 
term for the vorticity at TCA(II), the vorticity can be 
generated at this order. 



D. Generation of magnetic fields 



To evaluate the magnetic fields we substitute 6v 



{-tb)i 

[i.e., Eqs. (59), (63), (65), and (67)] into the source term 
for magnetic fields [i.e., Eq. (40)]. It is easy to see that 
the source term at TCA(I) is given by 



(2)i 



(72) 



where we used Eqs. (56), (71). This turns out to be pro- 
portional to the vorticity of the photon fluid. Indeed, we 
find 



+ 13 3e a(0) 



If one assumes a vanishing initial vorticity [see the ar- 
gument around Eq. (57)], magnetic fields cannot be pro- 
duced at TCA(I). 

The source term which we consider up to TCA(II) is 
expressed as 



iS],) = (5l5)(^) + (51j)(^^) 



(74) 



where a'--^'^\ and (^"^[^5)' are given by Eqs. (60), (68), 



and (59), respectively. Note in passing that ('^'y^^J,)'^)' is 



computed as 



Hence, we get 

2a^n 



W(3 + 2i?(o)) 



{1,1) 

1 + i?(o) "^(76)* 



+- 



a 



(0) 



(75) 



(0) 



1 + 13 3e 

(2)i I ijk 



1 + ^(0)°^? ^6 ^^(lb)k 



.(76) 



where we used Eq.(68). Intrinsically second-order per- 
turbations appear only as w^^^* in the source term, which 
is governed by Eq. (70). We can check that Eq. (76) is 
the same as the result obtained in Ref. [11]. Note that 
the vorticity is the general rclativistic effect because it 
vanishes in the limit of 7i ^ 0. Taking the Newton limit 
(H ^ 0), only the last term in the right-hand side of Eq. 
(76) remains. At first glance, this result is not the same 
as one obtained in Ref. [14] where the tight coupling ap- 
proximation was employed in Newtonian theory. This is 
not surprising because photon's rest frame is not used as 
the background of the tight coupling approximation in 
Ref. [14]. Performing an appropriate coordinate trans- 
formation, therefore, we can show that our result agrees 
exactly with [14]. 

Now we are ready to give an estimate of the magnetic 
field generated from second-order cosmological perturba- 
tions at the recombination epoch. For clarity and sim- 
plicity, we focus on the magnetic fields on the horizon 
scales (A; ~ Ti). We find 7?. ~ (5 by using the first order 
Einstein equations and v ^ {Tl/k)S ^ Shy using Eq. (44) 

on the horizon scales, where 5 = 5p^^ /p^\ Finally, the 
magnetic field on the horizon scales is estimated to be 



10 



-29 



Gauss, (77) 



where we used a ~ 10~^, k ~ (100 Mpc)"^ ~ 0.64 x 
10-^° GeV, e ~ 0.3, nip ~ 0.94 GeV, or ~ 1.7 x 

103 (GeV)-2, and p^°^a^ = p^o - 2.1 x 10-5i(GeV)'*. 
This result is consistent with the result in Refs. [9, 11, 13]. 
According to [2] , this will be amplified enough to explain 
the present observed magnetic fields. 



IV. SUMMARY 

In this paper we have derived an analytic formula for 
the magnetic fields generated from second-order cosmo- 
logical perturbations in the pre-recombination era. Since 
first-order vector perturbations are not generated from 
inflation in the standard scenarios, we must go beyond 
the linear analysis to generate magnetic fields. Photons 
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and charged particles are strongly coupled via Thom- 
son scattering within the cosmic plasma, and hence the 
system behaves almost as a single fluid. In this single- 
fluid description magnetic flelds are never generated, and 
therefore the tiny deviation from the single-fluid descrip- 
tion is crucial here. Using the tight coupling approxima- 
tion (TCA) to treat the small difference between photons 
and charged particles, we have seen how magnetic flelds 
are generated from cosmic inhomogeneities. It was found 
that magnetic fields are not generated at first order in the 
TCA. Therefore, we conclude that magnetogenesis re- 
quires both the second-order cosmological perturbations 
and the second-order TCA. The resultant magnetic fields 
are expressed in terms of the vorticity and the product 
of the first-order perturbations. The latter can be com- 
puted by solving the linear Einstein equations, while the 
former is governed by the vorticity evolution equation, 
with the source term given by the product of the first 
order terms. Finally, we gave an order of magnitude es- 
timate for the magnetic fields generated on the horizon 
scales around the recombination epoch. The result was 
B ~ 10~^® Gauss. This magnetic flelds will be enhanced 
by the dynamo mechanism during and after galaxy for- 
mation to explain the present order of the magnitude of 
the observed magnetic fields [2] . Although our rough es- 
timate suggests that the current generation mechanism 
can work even on smaller scales, say Mpc, more careful 
investigation is needed to make clear whether it indeed 
works in a wide range of spatial scales. 

This work is the general relativistic version of [14] , in- 
cluding the effects of the cosmic expansion and metric 
perturbations. The qualitative conclusion derived in this 
paper confirms the basic results obtained from the fully 
nonlinear, covariant analysis of Ref. [11]. However, our 
formulation here is based on cosmological perturbation 
theory, which will be more useful for actual calculations, 
e.g., of the spectrum of the magnetic flelds. 

We have not included the effect of anisotropic stresses 
of photon fluids and the recombination process, though 
they will be equally important for magnetogenesis on 
small scales [12, 13, 15]. Taking into account these effects 
and computing the power spectrum of the generated mag- 
netic flelds require detailed numerical calculations, which 
are left for forthcoming publications. 



APPENDIX A: CONSISTENCY WITH REF. [11] 

We show that our result can be identifled to that ob- 
tained by using the covariant approach [11]. The flnal 
result in Ref. [11] is 



= — 6B''+a\B^ 



R 



1-/33 



e A{l + RY l-h/32 



DcP-y 
Pi 



e 2(1 i?) 1 /32 



(Al) 



Here A(i) is baryon's density fluctuation at TCA(I), 6 e 
Vau" is the volume expansion rate, a\ = h^hff'V (c^d) ' 
\Ohab is the shear, 



for a tensor, 

for a scalar function and 

Da ■■= h\Vb 



(A2) 



(A3) 



(A4) 



where Ua is photon's four velocity and hab = 9ab + UaUb is 
the projection to 3-dimensional hypersurface normal to 
u". 

Equation (Al) is the nonlinear result. For purpose of 
comparison we expand Eq.(Al) in terms of cosmologi- 
cal perturbations up to second order. We use the per- 
turbed metric (1), and photon's four velocity is deflned 
by Eqs.(4) and (5). Using our notation, the i component 
of Eq.(Al) is 



^(2)' + f 6|(0)s(2)i ^ 

3 



l-/?3 



e 4(H-i?("))2 I + /32 
e 2(l + .R(0))l + /32^(0) 



,(2)i 



n -(1) 
(7,1) DuPi 



(A5) 
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We neglect the shear which is higher-order in cosmologi- 
cal perturbations. First, we calculate the left-hand side. 
Using Eq.(A2) and 0^°^ = SH/a, we find 



3 



a 



= Voi3(')' + 2-B(2)^ 
^ " a 

= 1 ((S(2)i)'+3WB(2)^) 

= ^,{a'B('^^)'. (A6) 
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Second, wc calculate the second term in the right-hand 
side. Using Eq.(41) and Eq.(A2), we obtain 



(2)i 



e 2(l + i?(o)) I + /32 pio) 



2(l+i?(o))l + /32 ^(c 



(0) 



,(2)i 



mp -R^") 1 - /j^ 
e 1 + ^(0) 1 + /32 



2n 



-co 



(2)i 



(A7) 



Third, we calculate the first term in the right-hand side. 
Using Eq.(A4) we get 



l-/?3 



(7,1) D^p^ 



(1) 



e 4(l + i?(0))2 I + /32 

'mp m 1 - /33 
. e (1 + E(0))2 I + /32 



4 ^ 



Finally, we substitute Eqs.(A6), (A7), and (A8) to 
Eq.(A5), and then we get 



e (l + i?(0))2 1 + /J2 



xae 



+ 



1 



ae 1 + i?(o) 1 + /32 
Using Eq.(60), Eq.(A9) can be recast in 



1+/3 3e 



■a p-y 



■'(76)/= 



. (AID) 



. (A8) This result is the same as Eq.(76). 
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